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Abstract 

We solve the Vlasov equation for the longitudinal distribution function and find 
stationary wave patterns when the distribution in the energy error is Maxwellian. 
In the long wavelength limit a stability criterion for linear waves has been ob- 
tained and a Korteweg-de Vries- Burgers equation for the relevant hydrodynamic 
quantities has been derived. 
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1 Introduction. 



Nonlinear wave interaction in high energy synchrotrons has recently received a great deal of 
attention (see e.g. [|1|, 0, 0), since it has proven its importance for understanding a variety 
of phenomena in high intensity beams. 

Perhaps, the simplest problem to study is the evolution in longitudinal direction only of 
a intense coasting beam influenced by a broad-band resonator type impedance. This model 
exhibits a surprisingly vast variety of interesting features, part of which have already been 
experimentally observed and theoretically investigated [Q, 0, 0. Different types of beam 
equilibria can be detected due to the collective (nonlinear) interaction between beam particles 
and resonator waves, the latter being induced by the beam itself. Solutions describing 
similar types of plasma equilibria [Bernstein- Greene-Kruskal (BGK) modes] are well-known 
in plasma physics |^ . Structures of arbitrary shape can be formed in the nonlinear stationary 
regime, which substantially depend on the type of the initial velocity distribution. 

It is the purpose of the present paper to apply techniques borrowed from plasma physics 
to study nonlinear patterns in coasting beams that are in close analogy with BGK modes. 
In sections 3 and 4 we solve the Vlasov equation by expanding the distribution function in 
a power series of the resonator potential 0, and in the case of initial Maxwellian energy 
error distribution we obtain an equation, describing the evolution of stationary waves on the 
resonator. In section 5 we find a stability criterion for linear waves in the long wavelength 
limit and derive a Korteweg-de Vries-Burgers equation for the beam density, current velocity 
and resonator voltage. 



2 Model. 

We consider the longitudinal dynamics of a high energy stored beam governed by the set of 
equations @, 0: 

dV di , , 

I(e;T)= ldvvf{e,v;T)^ (2.3) 

The first equation (2.1) is the Vlasov equation for the longitudinal distribution function 
f{9, v; T) of an unbunched beam, while the second equation (2.2) governs the variation per 
turn of the voltage V{9] T) on a resonator. All dependent and independent variables, as well 
as free parameters in equations (2.1-3) are dimensionless and have been rescaled according 
to the relations: 

1 = cogt ; v = = 1 -\ ; CO = — , (2.4a) 

uj^ dt uj^ Wo 
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ZC^ TT 

Here a;^ is the angular revolution frequency of the synchronous particle, AE is the energy 
error, ojr is the resonant frequency, Q is the quality factor of the resonator, TZ is the resonator 
shunt impedance and po is the uniform beam density distribution in the thermodynamic limit. 
Furthermore 

.0^-^ (2.5) 

is the proportionality constant between the frequency deviation of a non synchronous particle 
with respect to the synchronous one, while rj = um — 77^ (c^m - momentum compaction 
factor) is the phase slip coefficient. The voltage variation per turn V{6\ T), the beam current 
I{0;T) and the longitudinal distribution function f{9,v;T) entering equations (2.1-3) have 
been rescaled as well from their actual values Va{9;T), Ia{6;T) and fa{9,v;T) as follows: 

Va = 2euJspQ--fnV ; /„ = etOsPoI ; fa = Po/- (2.6) 
Prom the Vlasov equation (2.1) it is straightforward to obtain the continuity equation: 



_9_ /■ , „ d 
df 

which will be needed for the exposition in the next section. 



j dvf+—j dvvf^Q, (2.7) 



3 Solution of the Vlasov Equation. 

Let us now try to solve the Vlasov equation by the simple separation of variables ansatz: 

f{e,v-T)=g{v)^lj{e-T). (3.1) 
Substitution of (3.1) into the continuity equation (2.7) yields: 

i+oi = 0, (3.2) 

where 

_ !dvvg{v) 

Sdvg{v)- ^^-^^ 
The Vlasov equation (2.1) with (3.1-3) in hand can be further transformed to 

djj^ _ XVi/j dg 
'89 ~ g(n-v)du' 

The separation of variables ansatz (3.1) implies 
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dg Vt — V 
dv al ^' 

which leads to the well-known equilibrium Maxwell-Boltzmann distribution: 

\2- 



ij{e;T) = Zexp 



exp 



at 



v{e-T) 



de ' 



(3.4a) 

(3.4) 
(3.5) 



where 



27r 



Z 



-1 



d6 exp 



at, 



(3.6) 



The solution (3.4-6) suggests further generalization [Q of the separation of variables 
ansatz (3.1) 



f{e,v-T) = Y.g,{v)^\e-T). 

k=0 

Instead of equations (3.2) and (3.3) we now have 



where 



n{6;T) 



(3.7) 



(3.2a) 



k=l 



dvvgk{v) 



dvgk{v) 



(3.3a) 
(3.8) 

n{9; T) = const, (3.9) 

which will be proved a posteriori to hold and substitute (3.7) into the Vlasov equation (2.1). 
Taking into account (3.2a) we obtain: 



In order to determine the yet unknown functions gk{v) we make the assumption: 



{v - n)f:kgk{v)^^'\e- T) + Af:%M/'(^; T) = 0. 

k=l k=o ^"^ 

Equating coefficients in front of powers of Lp yields the following recurrence relation 



(3.10) 



{v-n){k + l)gk+i{v) = -A 



dgk{v) 
dv 



k + l 



or 

9k+i{v) = 

where we have introduced the operator 

V = 

Vt — V dv 

Noting that the formal solution of the recurrence relation (3.11) has the form 



d 



9k[v) 



k\ 



we finally arrive at the general solution of the Vlasov equation 



k=0 



What remains now is to verify the condition (3.9). It suffices to note that 



A* 



A;! 



Ak = YT dvvV^go{v) = — dv 



k\ 



dv 



V d 
Q — V dv 



k\ J (n-vY 



and similarly 



Thus 



dv 



-D'-'goiv). 



k\ J {n 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



Ak = nBk, (3.15) 

which proves equation (3.9). 

Clearly the solution (3.14) is uniquely determined by the generic function go{v). The 
simplest choice is when go{v) is the Maxwellian (3.4), that is go{v) itself is an eigenfunction 
of the operator T> with an eigenvalue a"^ [c.f. equation (3.4a)]. In this case we immediately 
recover the distribution (3.1) with (3.4-6). 



4 Nonlinear Stationary Waves. 



In order to derive an equation for the potential (p{9; T) we insert (3.1) and (3.4-6) into (2.2) 
and obtain: 



h h UJ 



ZQ 



d_ 



exp ( — 



(4.1) 
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Making use of relation (3.2a) we cast equation (4.1) into the form 



d^(f d'^ip r,d(p d 
dT^ ^dT^ dT dT 



exp 



Integrating once equation (4.2) with due account of the initial condition 



we obtain 



^ + 27^ + a;V = 2^^' 1-exp 



(4.2) 



(4.3) 



(4.4) 



Expanding the factor in square brackets on the right-hand-side of equation (4.4) around the 
stationary solution <^s = yields 



^ + 27^ + 0;^ = 



-(^ 1 + 



+ 



2(72 Qat 



+ ...]. 



(4.5) 



Above the transition energy 7^ > 7^ {^t = ctjw the parameter A is negative, so that 
two cases can be distinguished. Defining 



\X\ZVt^ 



UJq — UJ 



(4.6) 



we can state the two cases mentioned above in a more explicit way: 

Case I: Provided ujq > 0, equation (4.5) can be transformed to a damped Duffing equation 
with an additional quadratic nonlinearity 



dT 



2at 



1^1 ^' 



Case II: For a;o < equation (4.5) takes the form 



\uJo\(p 



>?Z^ 



1^1 / 



3(7,2 



(4.7) 



(4.8) 



In the limit 7 — > equation (4.8) can be solved when neglecting the cubic term. The result 
is: 



3|a;o|cr; 



X^ZQ? cosh^ 



(4.9) 



This is a drifting hump- like structure that is well-known as a solitary wave of the Korteweg-de 
Vries (KdV) type. 
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5 The Korteweg-de Vries-Burgers Equation. 

The exact solution of the Vlasov equation obtained in the preceding sections was found based 
on the stationary wave condition given by the continuity equation (3.2). In order to provide 
a more general treatment of the problem we introduce the new coordinates and variables 
along with the moving beam particles 

z = 9-T ; u = v-l. (5.1) 
Then the basic equations (2.1-3) can be written as: 

df . df df 
rP'V BV rl r 

-d^-^^^^ ^-d-zJ ^^^^ + ^^-^^ 
Let us now pass to the hydrodynamic description of the longitudinal beam motion. The gas 
dynamic equations read as 

■g^ + g-jm-O, (5.4) 

m+V^^XV- (6.6) 

d^V ^ dV dF dF 

__2,_ + .V=g^-^. (5.6) 

where 

F{z;9) ^ J duf{z,u;9) ; F{z;0)U{z;0) ^ j duuf{z,u;9). (5.7) 
Obviously the stationary solution of the gas dynamic equations (5.4-6) is given by 

Fo = l ; Uo^O ; 1^0 = 0. 
The dispersion law of linear waves of the form 

(F, U, V) - {Fl, Ul, Vl) exp [^ {Q9 - kz)] 
is governed by the following equation 

where Z{k) is the well-known impedance function 
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In the long wavelength limit (small k) the dispersion equation (5.8) has two roots given by 
the expression 

^1,2 = ^ ± v^A2 + 4Aa;2 + 4a;V2^ , (5.10) 

which are real below transition energy. However, the situation when the energy of the 
synchronous particle is above transition energy is different. The solutions (5.10) to the 
dispersion equation are real, provided 

|A|<Ai ; |A|>A2 ; Xl,2 = 2iu^(^lT^/l-(r^y (5.11) 
An instability occurs when Qi^2 are complex, that is when 

Ai < |A| < Aa. (5.12) 

In what follows we will study the case when our system is linearly stable, that is either below 
transition energy or in the stability region (5.11). 

The solution of the dispersion equation in the long wavelength limit suggests that new 
scaled coordinates should be introduced 0, |^ 

a = ^e{z -a9) ; x = e^^^^, (5.13) 
where e is a formal small parameter. Then the gas dynamic equations can be rewritten as 

dF d dF 

du ,au au aidF 

(5,6) 



da"^ da da dx 



where 



V = ^eV ; 7 = v^7o, (5.17) 

cjV-Aa-A-cjV^ = ; (^a = (5.18) 
Assuming the perturbation expansions: 

oo oo oo 

F=l+^e™F^ ; f/=^e-[/^ ; V = Y.^^'Vm (5.19) 

m=l m=l m=l 

for the first and second-order terms in e we obtain respectively 
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oa oa afj — 

or 

/7i(cr,x) = «Fi(a,x)+G(x), (5.21) 
where G{x) is a generic function of the variable x, and 

Ehminating F2, [/2 and V2 from equations (5.22) we finally arrive at the Korteweg-de Vries- 
Burgers equation 

(5.3) 



dx da da^ da"^ dx 



where 



^-^^ ; h^Y^r (5-246) 

Zauj'^ — A A — lauo'^ 

It is important to note that a~^Ui and A((T^ — a^)~^/ daVi satisfy exactly the same equation 
(5.23). 

Similar Korteweg-de Vries-Burgers equation in the case below transition energy has been 
recently derived by A. Aceves employing the method of multiple scales 0. 



6 Concluding Remarks. 

We have studied the longitudinal dynamics of a high energy coasting beam moving in a 
resonator. The coupled Vlasov equation for the longitudinal distribution function and the 
equation for the resonator voltage have been solved by closely following the method of 
Karimov and Lewis |^. The key point of this method consists in the representation of the 
distribution function as a power series in the resonator potential. Further self-consistent 
stationary wave patterns have been found in the simplest equilibrium case of Maxwellian 
distribution in the energy error. 
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In the long wavelength (small wavenumber) limit a stability criterion for linear waves 
has been obtained and a Korteweg-de Vries-Burgers equation for the relevant hydrodynamic 
quantities has been derived. 

An important (and interesting) extension of the results obtained here involves the longi- 
tudinal dynamics of a bunched beam. These will be reported elsewhere. 
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